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Abstract 

The shell structure of magic nuclei far from stability is discussed in terms of 

the self-consistent spherical Hartree-Fock-Bogoliubov theory. In particular, the 

^ I sensitivity of the shell-gap sizes and the two-neutron separation energies to the 



choice of particle-hole and particle-particle components of the effective interaction 
is investigated. 

1 Introduction 

Shell effects in nuclei lead to characteristic irregularities in otherwise smooth behaviour 
of global nuclear observables when they are analyzed as functions of Z, N, and A. This 
applies in particular to nuclear masses, or nuclear binding energies, which conspicu- 
ously show that the closed-shell nuclei are bound stronger than open-shell systems [1-3] . 
Rather simple one-body potentials can account for the values of particle numbers where 
strong shell effects occur [magic numbers, equal to 8, 20, 28, 40, 50, 82, and (for neutrons 
only) 126]. 

Through the binding energies, many other nuclear structure properties depend on 
the strength of shell effects (for discussion of various binding-energy relations, see Ref. 



[H, Sec. 2). This is particularly true for the particle-emission lifetimes and/or particle- 
capture probabilities which dramatically depend on accessible energies. In stable nuclei 
these dependencies are well studied and analyzed. However, for nuclei far from stability, 
the relevant experimental data are scarce or simply nonexistant. For example, the 
nucleosynthesis of heavy elements in supernova explosions depends crucially on the 
values of beta decay rates and the neutron capture cross sections in very neutron-rich 
nuclei which do not exist in terrestrial conditions [§, |^. A prediction of nuclear shell 
effects far from stability, based on the experience gained by studying stable nuclei, is an 
important and difficult challenge to nuclear structure theory. 

In the present contribution, we aim at quantitatively describing the magnitude of 
shell effects in nuclei close to the neutron drip line. The description is based on the 
self-consistent Hartree-Fock-Bogoliubov (HFB) theory with effective interactions fitted 
to known nuclei. In the particle-hole (p-h) channel, we employ the Skyrme interaction. 
In the particle-particle (p-p, pairing) channel, several contact forces are used. 

All details of the method follow Ref. 0, where the HFB theory in spatial coordinates 
has been introduced, and the Skyrme parametrization SkP has been obtained by a fit 
to properties of several magic nuclei and to the chain of the tin isotopes. In the present 
study, we also use the HFB method with other standard Skyrme-force parametrizations 
such as SHI |[ and SkM* §. 

In drip-line nuclei, special attention should be paid to the description of pairing 



correlations. In particular, the condition |T0|, W^: 



Sn ^ -A„ - A = 0, (1) 

which requires that at the drip line the neutron separation energy, Sn, vanishes, shows 
that pairing effects are crucial for determining the position of the one-neutron drip 
line. From Eq. (|l]) it follows that properties of the mean field, characterized by the 
Fermi energy. A, (determined by the p-h component of the effective interaction) and 
those of the pairing field. A, (determined by the p-p part of the effective interaction) 
are equally important far from the beta-stability line. Consequently, the choice of the 
pairing interaction is vital for analyzing shell effects in drip-line nuclei. 

In the present study, results obtained with several p-h and pairing interactions (see 
Section |^) are confronted with the data. A comparison of theoretical predictions for 
the two-neutron separation energies, S2n, with experimental values is given in Section ||, 
while in Section ^ the sizes of calculated neutron shell gaps for very neutron-rich nuclei 
are discussed. 



2 Pairing interaction 

2.1 Magnitude of pairing correlations 

Odd-even staggering of nuclear masses, i.e., the fact that nuchdes with odd numbers of 
protons or neutrons are found to have smaller binding energy than the arithmetic mean 
of the binding energies of the two neighboring even-even nuclei, is a direct consequence 
of nuclear pairing correlations. The strength of pairing interaction is often adjusted so 
as to reproduce the odd-even mass difference 

5n = ]^{B{N - 1, Z) + B{N + 1, Z)) - B{N, Z). (2) 

In particular, in Ref . , a certain linear combination of the parameters of the Skyrme 
force has been adjusted to fit the values of the odd-even mass difference in the tin 
isotopes. The resulting Skyrme interaction, SkP, describes both p-h and p-p channels. 
In Table |l] we present the neutron odd-even mass difference for proton-magic nuclei 
obtained in the HFB+SkP approach. For every magic proton number, calculations were 
performed for one odd isotope near the bottom of the stability valley and for its two even- 
even neighbours. In an odd isotope, a particle state was blocked (by the method given 
in Ref. |0]), and then the lowest energy was selected among different configurations. The 
experimental values of S^^ , determined from experimental binding energies of Ref. fT^ , 



are also given in Table |T]. In Table |^ we present the proton odd-even mass differences 
with all definitions and methods mutatis mutandis analogous to those used for neutrons. 
In Tables |l| and ^ we also present the values of experimental and calculated pairing 
gaps. The former are derived from finite-difference formulas involving experimental 
binding energies |T^ , and the latter are the average values defined according to Ref. |0] . 



As discussed in Ref. [|1], the experimental pairing gaps are, in general, larger than the 
odd-even mass staggering (0). 

Tables ^ and ^ show that the theoretical values of 5^ and 5-p are generally much 
smaller than the experimental values. The only exception is the odd-even mass difference 
for ^^^Sn which was actually used when fitting the SkP force. Based on this result, one 
may be tempted to refit the force parameters to increase the average magnitude of 
pairing correlations and improve the agreement between experimental and theoretical 
values of 5. However, a comparison of pairing-type two-body matrix elements of the 
SkP force with experiment ally- deduced pairing matrix elements []T^ indicates that they 



have, on the average, correct magnitude. We are therefore led to the conclusion that it 
is the blocking mechanism that is responsible for too low values of 5-^. 



In fact, the strength of pairing interaction in nuclei is never far away from the 
critical value at which, in the mean field approximation, the static pairing correlations 



vanish |T5[ . This situation precludes the use of the blocking approximation, because the 
critical value of pairing strength is rather large in many-quasiparticle configurations. For 
odd nuclei, particle number projection should, in principle, be used to describe pairing 
correlations near the phase transition region |jl6| . 



On the other hand, a comparison of experimental and theoretical (average) pairing 
gaps is less conclusive than a direct comparison of nuclear masses entering the odd- 
even staggering (|). This is so, because the staggering should specifically depend on 
the last occupied orbit near the Fermi energy, while the average pairing gap is a global 
property depending on all single-particle states. Until a better theory of pairing corre- 
lations in odd nuclei is applied, we do not have the correct tools which would allow the 
incorporation of experimental information on the odd-even staggering or on the experi- 
mental pairing gaps in a detailed adjustment of forces in the p-p channel. Therefore, in 
the present study, we consider even-mass nuclei only, and consequently, we use average 
pairing gaps. An and Ap, as qualitative measures of the pairing correlations. 

The last columns of Tables |l| and ^ display proton and neutron average pairing gaps 
at proton and neutron magic numbers, respectively. Neither iV=40 nor Z=40 shell gaps 
in this calculations are large enough to destroy pairing correlations. Static pairing is also 
obtained for the A^=20 and 28 systems. This suggests that the magnitude of pairing 
correlations predicted by SkP is slightly too large in light nuclei. However, this is a 
rather subtle discrepancy, because at particle numbers 20 and 28 (both for protons and 
neutrons) the total energies of paired and unpaired solutions are very close. For the 
protons, the unpaired solution is slightly lower in energy while for the neutrons, the 
paired solution is favoured. 

2.2 Pairing interactions 

In the present paper we use three kinds of pairing forces. Firstly, we present the HFB 
results with the Skyrme force SkP. As discussed in Refs. [7,17-19] and in the previous 
Section, the Skyrme parametrization can be adjusted in such a way that for states near 
the Fermi energy it has attractive pairing matrix elements. The SkP interaction fulfills 
this requirement. Secondly, we employ a simple contact interaction [20-22]: 

V'{r,r') = Vo6{r-r'), (3) 



and, thirdly, the contact interaction depending on the total (isoscalar) density of parti- 
cles p(r) [23-26], 

V'p{r, r') = {V, + V^p^) 6{r - r'). (4) 

Since the Skyrme force is also a contact interaction, the use of the SkP parametrization 
in the p-p channel amounts to adding momentum- dependent and spin-exchange terms 
to simple forces given in Eqs. @ and (^). 

The form of the density dependence used in Eq. (H) is completely analogous to the 
standard Skyrme force, and hence we use an analogous notation. Parameters Vq, V3, 
and 7 should be chosen in such a way that the interaction is strongly attractive outside 
the nucleus (i.e., Vq should be large and negative) and weak inside. In this way, the 
pair density becomes concentrated at the surface region while for the pure contact force, 
Eq. (Q), it is spread throughout the nuclear interior [24-28]. 

Parameters of phenomenological pairing interactions depend on the features of the 
mean field in the p-h channel, and especially on the density of single-particle states near 
the Fermi level. Therefore, for every parametrization of the Skyrme force used in the 
p-h channel, one should use different values of parameters defining pairing interaction. 
In the present work, we have adjusted the strength, Vq, of the contact force (0) in the 
following way: For every version of the Skyrme force used in the p-h channel, spherical 
HFB calculations (without blocking) were performed for ^^^Sn (see Section 2TT). Then 



the values of Vq were fixed by requiring that in each case the average pairing gap. An, 
was equal to 1.256 MeV which is the value obtained in the HFB-(-SkP calculations (see 
Table ^. When the contact interaction (|^) is used in the p-p channel and the SIII, 
SkM*, and SkP parametrizations are used in the p-h channel, we denote the resulting 
combined forces by SIII'', SkM'^, and SkP'^, respectively. 

Since the density dependence in Eq. (|) introduces two additional parameters, for 
the purpose of the present study, we have decided to use the same values of Vq and 7 as 
those used in the p-h channel, and then to fit V3 so as to reproduce the value of An in 
-^-^^Sn. These forces are denoted by SIII'^'', SkM'''', and SkP'''', respectively. Parameters 
of pairing forces used in the present study are listed in Table ^. In all cases, the phase 
space used, when solving the HFB equations, was chosen according to the prescription 
given in Ref. 0]. 

Average pairing gaps for semi-magic nuclei listed in Tables |I] and ^ are plotted in 
Figs. ID and |^ for neutrons and protons, respectively. Open and full symbols refer to 
pairing forces (|^) and (H), respectively. The results obtained with the SkP force agree 
quite satisfactorily with experimental values of odd-even staggering, apart from the 
neutrons in |?Ni28 where all the forces considered fail. Neutron gaps are, in general. 



reproduced better than the proton gaps which for some forces are often much too small. 
Therefore, for protons one sometimes has to use a larger interaction strength Vq than 



for neutrons p2|, ^. Here we are mostly interested in neutron properties and therefore 
we use a common value for both types of particles. Moreover, the SkP parametrization 
does not make such a distinction and, nevertheless, it provides a correct description of 
neutron and proton pairing properties simultaneously. 

Density-dependent pairing interactions work better in light nuclei. On the other 
hand, the pure contact force seems to work better for heavier systems. No general 
pattern appears as far as differences between various forces used in the p-h channel 
are concerned. Most likely, these differences are related to details of the placement of 
single-particle levels near the Fermi energy which significantly vary from force to force. 

As expressed by Eq. (P, to obtain a reasonable description of drip-line nuclei, pair- 
ing correlations have to be considered. However, compared with many detailed studies 
of Skyrme forces in the p-h channel, much less effort has been devoted to adjusting 
properties of pairing interactions. In stable nuclei, pairing correlations can be incorpo- 
rated to some extent by using the BCS approximation, but when departing from the 
beta-stability valley, use of the HFB method is a must |^, 2^. A simple one-parameter 



fit, carried out in the present study, makes it possible to reproduce an overall magnitude 
of pairing correlations in the few cases discussed. The forces listed in Table ^ should 
be viewed as zero-order approximations, enabling us to perform the HFB analysis with 
standard Skyrme interactions. Certainly, more extensive fits of contact forces (possibly 
with additional terms included) are necessary to obtain better overall agreement with 
experiment. In this respect, the relative success of the pairing SkP force may serve as 
an indication that such an improvement is indeed possible. 

3 Two-neutron separation energies 



Experimental two-neutron separation energies [0, S2n=B{N, Z)-B{N — 2, Z), are plot 



as 



ted in Fig. |^(a) for five chains of proton-magic isotopes. Values quoted in Ref. [^ 
"from systematics" are denoted by open symbols. Below the S'2„— 10 MeV line, experi- 
mental data are not available, and dramatic extrapolations in A^ are required to reach 
the region of the r-process path (S'2n~2 MeV) or the two-neutron drip line (5'2„=0). 

Among proton-magic, medium-heavy and heavy nuclei, there exists information 
about six closed neutron shells. Namely, the experimental data are available for two 
calcium isotopes at N=20 and A^=28, two nickel isotopes at N=28 and A^=40, and the 
zirconium and lead isotopes with A^=50 and A^=126, respectively. In all these cases. 



the values of S2n exhibit sudden jumps of varying magnitude [0. The size of the shell 



gap can be directly related to a difference between two-neutron separation energies, 

Si - Sl^^ = -B{N -2,Z) + 2B{N, Z) - B{N + 2,Z), (5) 

i.e., is related to the second derivative of the binding energy. Hence, in order to establish 
the size of the neutron shell gap, the mass of the doubly magic nucleus and the masses 
of its two neighbouring even isotopes have to be known. Similar information can be 
obtained by analyzing quantities A2n introduced in [Q. However, they involve masses 
of four isotopes and, therefore, they are less convenient. It should be stressed that in 
our study we do not aim at subtracting the smooth background from the experimental 
masses, but rather use indicators (^ to compare theoretical and experimental masses. 
A very long chain of known tin isotopes starts and terminates around neutron magic 
numbers A^=50 and N=82. Unfortunately, so far, the masses of ^°°Sn and ^^^Sn have 
not been measured, and ^^Sn is probably proton-unstable. Therefore, for the purpose 
of the present study, we establish the size of the N=82 shell gap at Z=50 by taking the 
mass of ^^^Sn from the systematic trend proposed in Ref. [^ . 



The neutron subshell effect at A^=40 is rather weak, and the experimental S'2n-values 
for the nickel isotopes exhibit only a rather small discontinuity [see Fig. ^(a)]. Therefore, 
in the following, we concentrate on neutron shell gaps obtained from the nuclear masses 
around 2oCa, 28Ca, 28^1; 5oZr, ^gi^^; ^^^ 126-^^. The direct information on the isotopic 
dependence of the shell-gap sizes is even more restrained. Only a comparison between 
the shell-gap values for ^^Ni and ^^Ca suggests that the size of the N=28 shell gap 
decreases with increasing neutron excess. 

In Fig. |](b), we compare, with experimental data, the results of global mass cal- 
culations obtained in the finite-range droplet model (FRDM) of Ref. |3^. This is a 



microscopic-macroscopic model which is known to give a very good fit to experimental 
nuclear masses. (It is worth noting that the FRDM explicitly considers deformation 
effects.) Consequently, the two-neutron separation energies are also well reproduced in 
this model. Only in some specific regions of the periodic chart are some discrepancies 
present. For example, a sudden decrease in S2n around ^^Zr is not reproduced and also 
the trend predicted for light zirconium isotopes is slightly off. As far as the jumps of S2n 
at magic numbers, the FRDM reproduces the six data points fairly well. In particular, 
the best agreement is obtained for ^"^Zr; the largest discrepancy is seen in ^^Ca. 

At this point, one should note that any kind of correlation, absent in a magic nucleus 
but appearing in its two even neighbors, would lead to a decreased size of the shell gap 
(H). Indeed, in such a situation, the binding energy of the magic nucleus is expected 



to be well described by an uncorrelated state, while the correlations in the neighboring 
nuclei would lead to a reduced value of >S'|^— S"^"*"^. Consequently, if correlations (e.g., 
fluctuations due to pairing vibrations, shape oscillations, etc.) are missing in the model 
description, the differences S"^— S'^^^ can easily be overestimated. Of course, this kind 
of reasoning tacitly assumes that the effects on shell-gap sizes originating from the pure 
mean field are correctly described which certainly need not to be valid. 

Predictions of the FRDM concerning the isotopic dependence of shell-gap sizes are 
very clear. With increasing neutron excess, the gap at A^=28 decreases (too much as 
compared with experiment), while the gaps at A^=50 and 82 stay practically the same. 
In practical terms, in the microscopic-macroscopic methods the shell-gap sizes are the 
direct consequences of simple parametrizations of the isotopic dependence of one-body 
average potential. Such dependences probably cannot be inferred, with sufficiently high 
precision, from fits to known nuclei, because these nuclei form isotopic chains which are 
not long enough. From this point of view, far-reaching extrapolations, based on two- 
body density-dependent effective interactions, have a better chance of success. Moreover, 
the results of microscopic- macroscopic methods are plagued with numerous uncertainties 
when they are applied to nuclei far from stability |^ . 



In Fig. §(c), we present results obtained by the spherical HFB-|-SkP calculations. 
Here, the overall agreement with data is similar to that obtained in the FRDM. Apart 
from slightly too high values of S2p obtained before the A^=126 shell closure in the lead 
isotopes, and small deviations for some zirconium isotopes (similar as in the FRDM), 
the isotopic dependence of the two-neutron separation energies is reproduced fairly well. 
The predicted size of the magic gap is slightly too large in ^''^Pb, slightly too small in 
^°Zr and ^^^Sn, and seriously underestimated in the light magic nuclei. Predictions 
concerning the isotopic dependence of the shell-gap sizes are markedly different than 
those of the FRDM. With increasing neutron excess, the shell gaps at A^=50 and 82 



decrease substantially. This effect has already been observed in Refs. ITT], ^ and is 
discussed in Section ||. 

After adjusting the strength of the contact force (^ to the odd-even mass staggering 
(Section |2.2|) , we are now able to perform the HFB calculations with other standard 
Skyrme parametrizations. The results for the SIII*^ and SkM*^ forces are presented in 
Figs, ^(a) and ^(b), respectively. They are compared with the results of the SkP*^ force 
shown in Fig. §(c). One should note that the results obtained with the SkP and SkP"^ 
forces are very similar, in spite of their very different properties in the pairing channel. 

As seen in Fig. ^, neither for SIII"^ nor for SkM'' is the agreement with data satisfac- 
tory. In particular, the shell-gap sizes at iV=50, 82, and 126 are strongly overestimated. 



and the values and the slopes of S2n are, in most cases, incorrect. At this point, one 
should note that a systematic few-MeV discrepancy in two-neutron separation energies 
leads to an integrated many-MeV discrepancy in total binding energies, cf. results for 
lead isotopes presented in Ref. [^]. 

The results shown in Fig. ^ illustrate a very strong dependence of the two-neutron 
separation energies on the force parameters. Although older forces, such as SIII and 
SkM*, can perform well in certain regions of Z and N , they do not give a satisfactory 
global reproduction of the data. On the other hand, a fairly good global agreement 
obtained with the SkP suggests that the improvement is possible, while a still better 



parametrization would be welcome (see Refs. [^, ^). Of course, forces which fail in 
reproducing the behavior with (A^ — Z) in known nuclei have little chance to perform 
better when going far from stability. For example, the almost constant values of 5*2,1 
obtained in the very heavy tin and lead isotopes with the SIII^ and SkM^ interactions, 
followed by strong shell effects at N=12Q and A^=184, do not seem very reliable. A 
detailed analysis of the force- dependence of results may give us valuable information 
on the relative importance of various force parameters. In the next Section, such an 
analysis is performed for the sizes of magic gaps in known and exotic nuclei. 

4 Magic shell gaps far from stability 

In the present Section, the Z-dependence of sizes of neutron shell gaps is analyzed. 
The shell-gap sizes are defined as in Eq. (^; six neutron magic numbers from A^=20 to 
126 are considered. In the Z-representation, neutron-rich nuclei appear at small proton 
numbers and the proton drip line is approached at large proton numbers. 

4.1 The role of deformation 

Experimental values of S^— S*^^^ as functions of Z are shown in Fig. ||(a). A con- 
spicuous feature present in the data is the appearance of well-defined maxima at magic 
numbers. This fact can be associated with the effect of deformation. Namely, the 
shell-gap sizes defined in Eq. (|) depend on the masses of even isotopes adjacent to a 
neutron-magic system. In nuclei which are not proton-magic, both neutron and proton 
shells are open. In such cases, deformation effects have to be considered. Therefore, the 
isotopic dependencies seen in Fig. |](a) illustrate not only the "pure" isotopic dependence 
of the shell effects but also the magnitude of deformation- like correlations. Of course, 
experimentally, one cannot easily disentangle the deformation and isotopic effects, and 
the complete theoretical description should take them into account simultaneously. On 



the other hand, extrapolations to nuclei far from stability are difficult enough, and one 
would wish to consider one problem at a time. Since, in the present study, we consider 
spherical shapes only, we must disregard the isotopic dependencies present in the data. 
Consequently, we concentrate mainly on the six cases of doubly magic nuclei discussed 
in Section ^. In Figs. ||(b) and ^(c), and also in the figures discussed below, these six 
values are denoted by full circles and, in addition, two points corresponding to the N=28 
systems are connected by a thick line. 

In spite of the fact that the FRDM results account for the deformation effects, 
this model does not reproduce the experimental differences 5'^— S'^'^^ in an entirely 
satisfactory way. The maxima at N=20 and N=Z=28 are not high enough and the 
maximum at ^^Ca is absent. The isotopic trend predicted for the light A^=50 isotones 
is correct, but the proton-number dependence in the heavy A^=50 isotones, and in the 
A^=82 and 126 isotones, is visibly too weak. One should, of course, bear in mind that 
we are now looking at very fine details of mass predictions. 

The HFB+SkP results presented in Fig. ^(c) show very clearly that sizes of spherical 
shell gaps strongly depend on the neutron excess. For all neutron magic numbers, 
the gap sizes decrease when approaching the neutron drip line; i.e., when the proton 
number decreases for a fixed value of A^. (Again, since these results are restricted 
to spherical shapes, one cannot expect to reproduce the experimental proton-number 
dependence seen in Fig. ^(a).) In the automated calculations we have performed for 
long chains of isotopes or isotones, in some cases, a poor convergence was obtained, 
due to crossings between single-particle or quasiparticle levels. Since, in this work, we 
discuss only general trends, we do not try to cure these cases. (This would require the 
full configuration blocking.) In Figs, ^^the converged points are connected by straight 
lines. 

4.2 Hartree-Fock results 

We are now in a position to discuss possible origins of the isospin dependence of the 
shell-gap sizes. First of all, by using the Hartree-Fock (HF) method in which the pairing 
correlations are neglected, we may study the role of the interactions in the p-h channel 
alone. In Fig. ^we present the HF results for three Skyrme forces, namely, SHI, SkM*, 
and SkP. 

Let us first discuss the force- independent, generic features of the predicted shell-gap 
sizes. For all the forces considered, the sizes of the N=20 to A^=82 shell gaps decrease 
with decreasing proton number. The effect is strongest for A^=28, and then weakens 
for A^=50 and A^=82. For A^=126, the gap size is predicted to be almost constant as a 
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function of Z, and even a reversed trend is obtained in the HF+SIII model. In all cases, 
a decrease of the N=28 shell gap between ^^Ni and ^^Ca is predicted. 

Systematically, the magnitude of the shell effect is too strong for the SIII and SkM* 
interactions. The sizes of the A^=50 and N=82 shell gaps are overestimated by as much 
as a factor of two. A part of this discrepancy can be attributed to missing correlations. 
As discussed above, they always decrease the shell-gap sizes; the influence of pairing 



correlations is discussed in Section ^^. In the HF+SkP model, the size of the N=28 
shell gap is too small. One may therefore anticipate that it will be much too small after 
the pairing correlations are taken into account. 

Overall magnitudes of shell-gap sizes can be clearly correlated with different values 
of the effective mass corresponding to the three studied Skyrme forces. Indeed, the 
single-particle level density g at the Fermi level can be approximated by [^, ^ 



where Nt is the particle number, kp is the Fermi momentum, and m* is the effective 
mass. For the SIII, SkM*, and SkP, the effective masses are ?7i*/r7i=0.76, 0.79, and 1.0, 
respectively. Therefore, the density of single-particle states increases when going from 
SIII to SkP, and hence the sizes of shell gaps decrease accordingly. 

Interesting conclusions can be drawn when comparing sizes of the A^=20 and N=28 
shell gaps. In experiment, the A^=20 gap is clearly larger (in ^°Ca) than the A^=28 gap 
(in ^'^Ca), which is the feature reproduced by the HF+SkP model. This means that 
the neutron if/2 orbital is closer to the upper major shell than to the lower one. On 
the other hand, the SIII and SkM* forces give similar sizes for both shell gaps, which 
may suggest that these forces place this orbital too low and have therefore too large a 



strength of the spin-orbit interaction (cf. discussion in Ref. [|34 



4.3 Hartree-Fock-Bogoliubov results 

By comparing results presented in Figs, ^(c) and ^(c), one may see that the pairing 
correlations dramatically influence the shell-gap sizes. In light nuclei, the gaps are 
reduced by about a factor of two, which gives substantial disagreement with experiment. 
One may attribute this effect to two factors: Firstly, the sizes of the A^=20 and N=28 
gaps obtained in the HF method are too small. This may suggest that in light nuclei 
the effective interaction should have an effective mass smaller than unity; this is the 
case for the SIII and SkM* forces. Since both predicted shell gaps at A^=20 and A^=28 
seem to be too small at the same time, the discrepancy with experiment should not 
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be attributed to an incorrect placement of the £7/2 shell, i.e., to an incorrect spin-orbit 
strength. Secondly, the pairing matrix elements of the SkP force are probably too strong 



in light nuclei. This has already been discussed in Section ^^1] in connection with the 
nonvanishing pairing correlations in the neutron-magic nuclei. Of course, the predicted 
strong pairing correlations can also result from reduced shell gaps, and both explanations 
strongly depend on one another. 

In heavier nuclei, pairing correlations bring the shell-gap sizes, predicted by the 
HF-l-SkP model, very close to experimental data. They also induce a dramatic quenching 
of shell effects in neutron-rich nuclei. For example, in the HF+SkP model, the N=82 
shell gap decreases from about 9 MeV to 5 MeV, with Z going down from 86 to 40. 
In the HFB+SkP method, the A^=82 gap decreases from 7 MeV to almost zero. As a 
result, there is almost no influence of the N=82 gap on the position of the spherical 
neutron drip line |]TT . 



Figures |^ and ^ display the HFB results for the shell-gap sizes obtained with the 
density-independent (||) and density- dependent (|) contact forces, respectively. The 
results are presented for the three Skyrme interactions discussed above. It can be 
seen that for the interactions SIII and SkM* which have low effective mass, pairing 
correlations are not sufficiently strong to bring the shell gaps in heavy nuclei down to 
the experimental values. In light nuclei, the sizes of the N=20 and N=28 shell gaps are 
again almost equal, and therefore either the former (for SIII^) or the latter (for SkM"^) 
values agree with experiment, but a consistent description cannot be obtained. 

In general, the shell-gap sizes obtained with the density-dependent pairing interac- 
tion are lower than those obtained with the pure contact force. The differences between 
different p-p interactions are significant. For example, the A^=82 shell gaps predicted 
in the SkM'' calculations are always greater than 5 MeV, while those obtained in the 
SkM'''' version exhibit a pronounced quenching when approaching the neutron drip line. 

5 Conclusions 

The results of self-consistent HF and HFB calculations with several Skyrme-type in- 
teractions demonstrate that the magnitude of shell effects at magic numbers strongly 
depends on force parameters. In the p-h channel, the main factor is the effective mass 
which determines the single-particle level density around the Fermi level. In the pair- 
ing channel, the density-dependence of a p-p interaction plays a very important role. 
Both effects combined, i.e., the increase in m*/m and the presence of density-dependent 
pairing, lead to the shell- quenching effects in nuclei far from the beta-stability line. 



12 



This result seems to be consistent with the recent analysis of solar-system r-process 



abundance distributions [0 



Our brief analysis indicates that a lot of systematic work still needs to be done to 
obtain a "universal" Skyrme-type interaction. For instance, the density dependence of 
the effective interaction invites a much deeper analysis. (In principle, the force should 



depend on both the isoscalar density, p„ + pp, and the isovector density, pn — Pp [^3 HI 
the latter might be important in the description of exotic nuclei with a significant neutron 
excess.) Also, the interplay between deformation and pairing in open-shell drip-line 
nuclei is still awaiting explanation. 
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Table 1: Experimental and calculated neutron odd-even mass difference S^^ and S-^ for 
selected proton-magic nuclei. Experimental pairing gaps A^^ (determined as in 
and calculated neutron and proton average pairing gaps An and Ap are also given. 



N 


Z 


rcxp 

On 


5n 


/\exp 


An 


Ap 


23 


20 


1.600 


1.301 


1.708 


1.584 


0.000 


31 


28 


1.194 


0.234 


1.446 


0.720 


0.000 


53 


40 


0.743 


0.448 


0.829 


0.712 


0.793 


69 


50 


1.312 


1.293 


1.378 


1.256 


0.000 


121 


82 


0.736 


0.634 


0.801 


0.737 


0.000 



Table 2: Same as in Table |l|, except for protons replacing neutrons and 



vzce versa. 



N 


Z 


5;^P 


Sp 


A^^P 


Ap 


An 


20 


17 


0.929 


0.496 


1.535 


1.191 


1.202 


28 


25 


1.147 


0.709 


1.540 


1.139 


0.605 


40 


31 


0.938 


0.303 


1.385 


0.900 


1.319 


50 


37 


0.997 


0.514 


1.356 


0.976 


0.000 


82 


59 


0.998 


0.664 


1.229 


0.931 


0.000 


126 


85 


0.661 


0.530 


0.846 


0.552 


0.000 
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Table 3: Parameters of pairing forces used in the present study. Columns denoted by S 
and Sp give parameters of forces (H) and (|[), respectively. In the latter case, parameters 
Vq and 7 are identical to those in the p-h channel. Different sets of parameters are fitted 
to different forces used in the p-h channel, as indicated in the first column. 



Force in the 


6 


6p 


p-h channel 


Vo (MeV fm^) 


Vo=to (MeV fm^) 


Vs (MeV fm^+T) 


7=a 


SIII 


-160 


-1128.75 


50715 


1 


SkM* 


-173 


-2645 


21591 


1/6 


SkP 


-185 


-2931.696 


24175 


1/6 
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Figure captions 



Figure 1: Average neutron pairing gaps, An, for proton- magic nuclei listed in Table |l|. 
Odd-even staggering parameter 6^^ is plotted as the experimental value. All forces were 
adjusted so as to give identical results for ^igSusQ. 



Figure 2: Average proton pairing gaps, Ap, for neutron- magic nuclei listed in Table 
Odd-even staggering parameter ^p""^ is plotted as the experimental value. 



Figure 3: Two-neutron separation energies, S2n, for proton-magic isotopes. In (a), 
experimental (full circles) and "from systematics" |]l2l (open circles) values are shown. 
In (b) and (c) compared with the experimental data are the results of the FRDM [^ 
and of the HFB-|-SkP model, respectively. 



Figure 4: Two-neutron separation energies, S2n, for proton- magic isotopes. The HFB 
results with the SIII*^ (a), SkM*^ (b), and SkP'^ (c) parametrizations (see Section pl^ are 
compared with the experimental data. 



Figure 5: Sizes of neutron shell gaps as functions of the number of protons. Experimental 
data (a) are compared with the results of the FRDM [^ (b) and of the HFB-|-SkP 
method (c). Deformation effects are included in the FRDM while only spherical shapes 
are considered in HFB-|-SkP. In (b) and (c), experimental data are shown only for magic 
nuclides and the thick line connects experimental points for two magic A^=28 systems. 



Figure 6: Sizes of neutron shell gaps obtained within the HF approximation for three 
parametrizations of the Skyrme force: SIII (a), SkM* (b), and SkP (c). Experimental 
data are shown only for magic nuclides and the thick line connects two experimental 
N=28 points. 
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Figure 7: Sizes of neutron shell gaps obtained within the HFB approximation for three 
parametrizations of the Skyrme force: SIII^ (a), SkM'' (b), and SkP'' (c). Contact 
pairing force (^ is used. The same experimental data as in Fig. ^ are shown. 



Figure 8: Sizes of neutron shell gaps obtained within the HFB approximation for three 
parametrizations of the Skyrme force: SIII'^'' (a), SkM^^ (b), and SkP'^'' (c). Density- 
dependent contact pairing force (^) is used. The same experimental data as in Fig. ^ 
are shown. 
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